Dynamical ultrametricity in the critical trap model 
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We show that the trap model at its critical temperature presents dynamical ultrametricity in the 
sense of Cugliandolo and Kurchan jl|] . We use the explicit analytic solution of this model to discuss 
several issues that arise in the context of mean-field glassy dynamics, such as the scaling form of 
the correlation function, and the finite time (or finite forcing) corrections to ultrametricity, that are 
found to decay only logarithmically with the associated time scale, as well as the fluctuation dissi- 
pation ratio. We also argue that in the multilevel trap model, the short time dynamics is dominated 
by the level which is at its critical temperature, so that dynamical ultrametricity should hold in the 
whole glassy temperature range. We revisit some experimental data on spin-glasses in light of these 
results. 
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PACS numbers: 75.10.Nr, 05.20.-Y, 02.50.-R 
I. INTRODUCTION 

Notable theoretical progress in our understanding of 
the ubiquitous aging phenomena in glassy systems was 
made possible by the recognition that the discussion of 
correlation and response functions requires two times: 
the waiting time t w and the total elapsed time t w + t. 
This appears very clearly in the framework of mean-field 
spin glass models, domain growth (coarsening models), 
or more phenomenological trap models Q . Although the 
basic phenomenology of all these models are rather sim- 
ilar, the underlying physical picture is completely dif- 
ferent. For example, aging in domain growth models is 
associated to the growth of a coherence length. In mean 
field or trap models, space is absent and aging is related 
to the structure of the energy landscape, but here again 
the intuition is completely different. In mean field mod- 
els, the system never reaches the bottom of an energy 
valley and there are no activated processes involved in 
the dynamics. Rather, the dynamics slows down because 
saddles with less and less 'descending' directions are vis- 
ited as the system ages. Conversely, in the trap model, 
activation is the basic ingredient of the model, and aging 
is associated to the fact that deeper and deeper valleys 
are reached as the system evolves. Dynamics in the lat- 
ter case in fundamentally intermittent: either nothing 
moves, or there is a jump between two traps. This must 
be contrasted with mean field dynamics which is continu- 
ous in time. However, mean field and trap dynamics can 
be shown to correspond to two successive time regimes 
in a particular class of models || . 

In spite of these important differences, many predic- 
tions are common to the latter two pictures, such as: 

• A short time singularity of the response function 
in the aging regime, which leads to an (aging) low 
frequency noise. 

• Non trivial violations of the Fluctuation Dissipa- 
tion Theorem (fdt), first pointed out within mean 



field models, but that also exists in trap models. 

• The possibility of rejuvenation and memory, which 
involves the existence of different degrees of free- 
dom with different time scales. 

Furthermore, a certain class of mean field models (that 
includes the Sherrington-Kirkpatrick model) has been 
conjectured to possess ultrametric dynamical properties, 
that very precisely reflect and encode the ultrametric na- 
ture of the static solution. This dynamical ultrametricity 
is associated to an infinite number of time scales (which 
diverges with the age of the system), in the following 
sense: if G(t2,ti) is the correlation function between 
times t\ < t 2 , then in the limit of large times: 

C(t 3 ,t 1 )=min(C(t 2 ,t 1 ),C(t 3 ,t 2 )), Vt 2 G [ti,t 3 ]. (1) 

This means that either t 2 is close enough to £3, and then 
no further dynamics takes place between t 2 and t 3 , or 
t 2 is close enough to t\ but then the age of the sys- 
tem hardly changes between t\ and t 2 . This property 
of the correlation function has been shown to hold for 
the aging solution of the dynamical (Mode Coupling) 
equations describing the dynamics of 'continuous' spin- 
glasses. This property is furthermore invariant under 
reparametrizations of time, where t — ► h(t) with an ar- 
bitrary monotonous function h(t). Testing whether or 
not dynamical ultrametricity also holds in realistic disor- 
dered systems is made difficult because this property is 
only expected in the limit of asymptotically large times, 
and corrections are expected on finite times. How large 
are these corrections? 

In this paper, we show that exact ultrametricity holds 
at the critical point of the single-level trap model (or ran- 
dom energy model) . We give the explicit form of the cor- 
relation function and discuss finite time (or finite forcing) 
corrections. Note that in this single-level trap model, the 
dynamics is ultrametric although the statics is not. The 
issue of finite time fdt plot is also adressed. We dis- 
cuss multi-level extensions of the trap model and argue 
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that dynamical ultrametricity should be generic at 'short 
times', i.e. at the beginning of the aging region. We 
show thermoremanent magnetization data that support 
this idea. The relation with 1// noise, already discussed 
in this context H], is recalled. 



II. THE MODEL 

The trap model, introduced in the context of aging in 
||,[| and further developed in Q , is one of the simplest 
soluble model exhibiting a dynamical glass transition. In 
this model, one considers a particle which is trapped in 
low energy states i of depth Ei > 0, where the Ei are ran- 
dom variables distributed according to p(E) = jre~ E / Tc . 
The dynamics is chosen to be activated: each particle 
stays in trap i an exponential random time, equal on av- 
erage to TQe Ei ' T . The quantity To is a microscopic time 
scale which we shall take as the time unit in the following. 
When the particle leaves the trap, it chooses at random 
a new one among all the others. As a consequence, at 
high temperature (T > T c ), the particle spends most of 
its time in the small traps, because the number of these 
traps (the entropic factor) dominates the Boltzmann fac- 
tor and the system equilibrates. On the contrary, for 
T <T Cl the Boltzmann factor is dominant and the parti- 
cle explores deeper and deeper traps, so that the system 
never equilibrates (in the limit of an infinite number of 
traps). In this regime, the dynamics ages: correlation 
and response functions are no longer time translation in- 
variant, but depend both on the waiting time t w and the 
total time t w + t. 



III. CORRELATION FUNCTION 

Correlation functions are useful tools to characterize 
the dynamics and to compare several models. The sim- 
plest (but nevertheless non trivial) correlation in the trap 
model is H(t w +t, t w ) defined as the probability to remain 
in the same trap during the time interval [t w , t w + 1] . As 
we are considering the infinite dimensional (or fully con- 
nected) model, this probability is also equal to the prob- 
ability P(t w +t, t w ) to be in the same trap at t w + t as at 
t w , since the probability to go back to the same trap van- 
ishes in the limit of an infinite number of traps. (This is 
not true in finite dimension. For instance, in one dimen- 
sion, H(t w + t, t w ) and P(t w + t,t w ) scale in a different 
way with t w ). 

The correlation function Tl(t w +t,t w ) has been calcu- 
lated in the general case in |Q. We shall now focus on 
the critical case T — T ri for which we have: 



for large t w . We note that this expression simplifies fur- 
ther in the limit t, t w — > oo. In particular, the exponential 
term vanishes since t w + 1 — u > t. Writing u = t w (\ — v), 
we get 



dv 



pog *„+ log (!-«)][« + 



TT ( 3 ) 



dz 



O 



logtw Jq z + 1 \(logt w ) 



(4) 



with the new integration variable z = This expres- 
sion is readily integrated to give: 



log(l + V) 
logt w 



(5) 



This relation shows that H(t w +t,t w ) is not a function of 
j-, at variance with the results that hold the whole low 
temperature phase T < T c 0. On the contrary, taking 
t ~ at^, we obtain, in the limit t w going to +oo: 



n(t„ +t,t w ) =C(a) 



with C(a) given by: 



C(a) = l-a (a < 1) 
= (a>l). 



(6) 



(7) 
(8) 



Note that C(a) is a monotonous decreasing function of 
a. The above result is only true in the limit t tq, as 
well as log (£u,/to) S> 1. Indeed, all the finite time results 
reported in this paper should be understood as first order 
terms in a 1/ log t w expansion. 

An important remark is that the correlation function 
n(i„, +t, t w ) is a function of a = log tj log t w that cannot 
be written as h(t w + t)/h(t w ). The latter ratio naturally 
appears (with an unknown function h) in the aging part 
of the solution of the dynamical equation corresponding 
to one step Replica Symmetry Breaking (rsb) mean field 
spin glass models. However, for full RSB models where 
dynamical ultrametricity is indeed expected, the correla- 
tion function is given by an infinite sum of contributions 
coming from different times sectors 0,^]: 



C(t w + t, t w ) — N ' C 



hj(t w + 1) 

hi (^u>) 



(9) 



where hi are unknown (monotonous) functions defining 
the ith time scale, and the Ci are monotonously decaying 
to zero for large arguments. 

A useful (but up to now unjustified theoretically) form 
for hi(t), that allows one to give some flesh to the above 
formula, is 0: 



du 



1 - exp [-(t w +t - u)j 
(logu)(t w + t-u) 



(2) 



hi(u) = exp 



/./,,■ 



< m < 1. 



(10) 
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It is easy to see that for this choice of hi, the time scale 
on which the ratio hi(t w + t)/hi(t w ) varies significantly 
is precisely t$. The choice \i = therefore corresponds 
to stationnary dynamics, whereas p, = 1 gives full aging. 
Now, take t = t° (with < a < 1) in Eq. (g) and take 
the limit t w — ► oo. All the sectors such that /ij < a 
have relaxed to zero, whereas the sectors corresponding 
to fj,i > a have not decayed at all. Introducing a contin- 
uum of different values of p, we find that the correlation 
function is given by: 

C{t w +tl,t w )= I dfip(p)C^l), (11) 

J a 

where p(p) is the 'density' of time sectors of order fc^ and 
C M (1) is the initial value of the correlation function in this 
sector. From this result, one sees that C(t w +t, t w ) indeed 
becomes a function of a = log i/ log in the long time 
limit. Therefore, interestingly, the superposition of an 
infinite number of subaging contributions defined by (|Io| ) 
naturally leads to a correlation function that depends 
on log tj log t w , for which the dynamical ultrametricity 
property is explicit. The critical trap behaviour corre- 
sponds to a uniform contribution of all time sectors, i.e. 
= 1, V/x. 



• If P < a, then one has, for large t\: 

t 3 -t 2 ~at% -btf ~at% (15) 

so that Ci = C(/3) and C* 2 = C 3 = C(a) < C{(3). 

• Assuming now j3 — a, we get C± = C3 = C{a), as 
well as C 2 > C 3 . 

As a result, we have shown that the relation C3 = 
min (Ci, C 2 ) always holds, which implies that dynamical 
ultrametricity is satisfied in this model. 

The appearance of dynamical ultrametricity can be 
considered as a signature of the existence of many time 
scales involved in the dynamics. This is indeed the case 
in the present model, even though there is no static ul- 
trametricity to account for this hierarchy of time scales. 
This property in fact arises naturally from the exact bal- 
ance, at the critical point, of the Boltzmann weight e E l T 
and of the entropic factor p{E): the probability for the 
particle to have a given energy (or equivalently, a given 
trapping time, in logarithmic scale) is essentially uniform 
on the interval [0, T c logt w ]. In other words, dynamical 
ultrametricity here is a consequence of the critical scale 
invariance. 



IV. DYNAMICAL ULTRAMETRICITY 

As recalled in the Introduction, Cugliandolo and Kur- 
chan have defined dynamical ultrametricity for the cor- 
relation function C if the following property is true: take 
three times t\ < t 2 < £3, then: 

C(t 3 ,t 1 ) = mia{C{t 2 ,t 1 ),C(t 3 ,t 2 )), Vt 2 G[ii,* 3 ]. (12) 

Let us show now that H(t w + t,t w ) at the critical tem- 
perature is ultrametric in the sense defined hereabove. It 
will be useful to introduce the following notations: 

n(t 2 ,ti) = Ci (13) 

U(t 3 ,t 2 ) = C 2 

n(t 3 ,h) = c 3 

From the monoticity of correlation functions, the inequal- 
ity C3 < min(Ci,C2) holds in general. We simply have 
to check that (at least) one of the two correlations C\ 
and C 2 is equal to C 3 . In order to take the infinite time 
limit, we need to specify how t 2 and £3 scale with t\. A 
natural parametrization is the following: 

t 3 -t x ~atl (14) 
We can now look at various cases: 



V. FINITE TIME ANALYSIS 

An interesting analysis was introduced by Cugliandolo 
and Kurchan in the context of the dynamical analysis 
of the Sherrington-Kirkpatrick (SK) model M. These 
authors made the assumption that, in the limit of large 
times, there exists a certain function f(x,y), not neces- 
sarily smooth, such that C3 = f{C\, C 2 ). We have shown 
in the previous section that such a function indeed exists 
in the present case and is given by f(x, y) = min (a;, y). 

A useful representation, proposed in [0, is to plot in 
the (Ci, C2)-plane the curves of constant C3 = f{C\,C 2 ), 
which reduces for our case to two straight lines (C\ = C 3 
and C 2 = C 3 ) at right angle. For finite times, the 
function / has to include a time scale as third argu- 
ment, so that C3 = F{C\,C 2 ;t\), where f(x,y) — 
lim tl _ too F(x, y; t\). The (Ci, C2)-plane representation 
is then a good way to visualize the convergence towards 
the asymptotic function f(C\,C 2 ). It has been used with 
numerical data to test if dynamic ultrametricity holds in 
realistic systems, with rather inconclusive results 
except in a few cases, like in a recent study of the 4- 
dimensional Edwards- Anderson model || . 

Let us apply this procedure to the critical trap model. 
In order to deal with finite time expressions, we shall 
come back to Eqn. H, restated as follows: 

log (1 + 5^7) 
C ^= log£ (16) 
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Inverting these relations so as to express t 2 and t% as 
functions of t\, C\ and C3, we can write an explicit ex- 
pression for F(Ci, C2; t\): 



C 3 = F(C u C 2 ;t 1 ) 



log(*r 01 +t- C2 (l-q 



-ClNl + C 2 



(17) 



In order to plot the constant C3 curves, it will be useful 
to express also C 2 as a certain function G(C\, C3; t{): 



io g (i-ir Cl )-iog(ir°' 3 -ti Ui ) 



-C 3 



C2 — G(Ci, C3; ti) 



log*! -log(l -t- Cl ) 



(18) 



The resulting plots in the (C±, C2)-plane are displayed in 
Fig. § for C 3 = 0.1,0.3,0.5,0.7 and for three different 
times (t\ = 10 5 , 10 10 and 10 15 ). Note that the conver- 
gence is very slow close to the infinite time singularity. 




F( C, C;tl) = C-^±<^p^l (19) 

logti 

For C > 0, this expression simplifies further at large 
times: 



F(C,C;h) ~C - 



log 2 



(20) 



Interestingly, the leading correction does not depend on 
C, and it is valid only if C is not too close to 0. Fig. ^ 
displays the plots of F(C, C; ii) — C, for the same values 
of ti as in Fig. [|. 




FIG. 2. Plot ofF(C,C;ti) - C versus C for the same val- 
ues of ti as in Fig. ^|. The dotted curve, corresponding to 
the infinite time limit f(C, C) = C is added for comparison. 
The departure from f(C, C) — C is almost independent of C, 
except for C close to 0. 



FIG. 1. Plot of constant C 3 = F(Ci,C 2 ;t 1 ) in the 
(Ci, C2)-plane. From left to right, curves correspond to 
C3 = 0.1,0.3,0.5,0.7. Each set of three curves shows 
the convergence with t\ towards the asymptotic function 
/(Ci,C* 2 ) = min(Ci,C 2 ): h = 10 5 (full line), 10 10 
(dot-dashed line), 10 15 (dashed line). 

Cugliandolo and Kurchan also introduced in JjJ the 
notion of correlation time scales through the representa- 
tion of f(C, C) versus C. As already mentionned before, 
f(C,C) < C in general. There may exist some special 
'fixed' points C* such that f(C*,C*) = C*. Each of 
these fixed points has been shown to be associated with a 
correlation time scale. If dynamical ultrametricity holds 
in a particular time sector, all C's belonging to a cer- 
tain interval [C',C"} are fixed points. In our case, ul- 
trametricity holds over the full correlation interval [0, 1]. 
But in our model, we can go beyond the infinite time 
analysis and quantify the convergence of F(C, C; ti) with 
t\ towards the asymptotic function f(C, C) = C . We 
find: 



This last result may also be interpreted in the frame- 
work of Fig. 0. For a given value of C3, the point 
C(t x ) defined by the relation F(C(t x ), C(ii);*i) = C 3 
converges to the infinite time right angle singularity 
d = C 2 = C* 3 as (see Eqn. poh: 



C(*i 



C 3 + 



log 2 
logti 



(21) 



This logarithmic correction may explain why it seems so 
difficult to observe the convergence towards the ultra- 
metric relation in experimental or numerical data, where 
only a few decades (usually between four and six) are 
available. 



VI. THE EFFECT OF 'SHEAR' 

The effect of an external 'shear' (or power injection) 
on aging was investigated in the context of mean field 
models in H, and in the context of the trap model in 
p0|. In both models, aging is interrupted by the shear 
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beyond a time scale r r which diverges as the shear rate 7 
tends to zero. In the model considered in ]To[ |, this time 
scale is given by: 



7 V 7 



(22) 



For r <C r r , the (power-law) distribution of trapping 
times is in a first approximation unaffected by the shear, 
whereas for longer times, the distribution decays expo- 
nentially. In the limit where the waiting time t w is much 
larger than r r , the dynamics of the model becomes sta- 
tionary and one finds, for t <C r r , the same result as 
above with t w replaced by r r : 



C(t + tyj, t W ) C 



log 17 



1 



logT r 



(23) 



As discussed in |8|,|9|], dynamical ultrametricity in this 
context manifests itself by the appearance of an infinity 
of time scales in the limit r r — > 00 (i.e. 7 — > 0): the time 
needed for the correlation to decay to a certain value c 
diverges as t}~ c (see Q for a further discussion). 



one sees a very slow convergence towards the asymptotic 
result for small values of C. 




C(t 2 ,t t ) 



FIG. 3. Plot of integrated response x(t2 , 1 1 ) versus C(t2,t\) 
parametrized by t2, for t\ — 10 3 , 10 5 and 10 10 from bottom 
to top. The local slope is equal to —X(t2,t\)/T c , and it con- 
verges very slowly towards the asymptotic value —1/T C for 
small C (see inset). T c is chosen as the temperature unit. 



VII. THE FLUCTUATION-DISSIPATION RATIO 

It is interesting to study the fluctuation dissipation ra- 
tio X in the trap model at the critical point. This ratio 
is defined as: 



TR(t2,h) 
dC(t 2> ti)/dti 



(24) 



where R(t2,ti) is the response of the system at time 
t 2 > fi to a small bias field applied at time t\ (see 
jj],[l2| for details). In the trap model, where the field 
only changes the trapping time of the starting site, one 
finds that the following relation holds in general: 



TR(t 2 ,h) 



dC(t 2 ,h) 
dt 2 ' 



(25) 



Using this result, one sees that in the 'liquid' phase 
T > T c where all two time functions only depend on time 
differences, X = 1: the usual Fluctuation-Dissipation 
Theorem (fdt) holds. In the glass phase, on the other 
hand, one finds that X(t 2 , t\) — t\jt 2 : the value of X is 
non trivial in the whole scaling regime t 2 — t\ ~ t\. Right 
at the critical point T = T c , one can express X as: 



X(t 2 ,tx) 



(*f - I) 2 



C = C(t 2 M) >0. (26) 



So for any fixed C > 0, X tends to 1 in the asymptotic 
limit t\ — * 00. We show in Fig. ^| the now famous plot 
of the integrated response versus C, which should yield a 
straight line of slope — 1/T C when the fdt holds. Again, 



VIII. THE MULTILEVEL TRAP MODEL AND 
DISCUSSION 

The simple trap model described above can be consid- 
ered as a one-step RSB model. In order to generalize the 
model to a full RSB one, it has been proposed in (and 
further studied in p3| ) to follow Parisi's procedure for 
the static solution of the SK model. Roughly speaking, 
it means that each trap is recursively subdivided into a 
new series of traps, in a hierarchical manner. Each level k 
of traps is characterized a certain overlap between states 
qk and by an exponential probability distribution of the 
energy barriers, with a critical temperature T c fc depend- 
ing on the level index k. The critical temperatures are 
related to Parisi's function Xk = x(qk) as T c fe = T/xk, 
and satisfy the relations T c fe < T* -1 . At any tempera- 
tures, the levels of the tree corresponding to q > qea(T), 
where qEA is the Edwards- Anderson parameter, are such 
that Xk > 1, so that these levels are equilibrated (i.e. 
T c fc < T). 

In the single level trap model, the correlation function 
H(t w +t,t w ) in the aging phase T <T C behaves at short 
times as: 



U(t w +t,t w ) ~ 1 



sin ttx ft 



7r(l — X) \t 



K t < tu 



(27) 



with x = ijr . In the multi- level model with a finite num- 
ber M of levels, the total correlation function is defined 
as: 
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M 

E 

fc=0 



g*: [n fc [t v 



Hk+l(t W + t, ty 



M 



= go + / J(gfc ~ gfc-l) life fa + t, (28) 
fc=l 

(ft being the fc" 1 level overlap, and Il^fa is the 

probability that the process has never jumped beyond the 
fc th layer of the tree between t w and t w + 1 with the con- 
vention that Ilofa +t, t w ) = 1 and n.M+ifa +t, t w ) = 
(see j7j for details). From this definition, we see that 
C(t w + t,t w ) is dominated at short times by the levels 
k with Xk close to 1, for which the short time singular- 
ity is strongest. (We assume that T < T c °, i.e. that at 
least one level is aging). Therefore, we expect to observe 
the dynamical ultrametricity associated to the level k* for 
which Xk* = 1 in the 'short' time regime logi/ logi„, < 1, 
before the t/t w regime associated to the levels k < k* 
sets in. (Note that if C(t w + t, t w ) is a function of t/t w , 
then the function f(x,y) defined above cannot be equal 
to min(x, y)). 

Interestingly, 'short time' dynamical ultrametricity ex- 
ists for the hierarchical tree model in the whole low tem- 
perature phase, but has no relation with the static ul- 
trametricity built in the tree structure which encodes 
Parisi's rsb solution. Thus the origin of dynamical ultra- 
metricity in the generalized trap model is again very dif- 
ferent from the dynamical ultrametricity found in mean 
field spin glass models, which in the latter case is deeply 
related to the Parisi function x(q) which encodes the 
structure of the tree. As mentioned in the introduction, 
the physical interpretation of aging in the two pictures 
are radically different, although some of the phenomenol- 
ogy is very similar. 

Aging experiments is spin-glasses have been inter- 
preted within the framework of the multi-level trap model 
in 0. The need for several levels comes from rejuve- 
nation and memory in temperature shift experiments, 
but also from the detailed shape of the thermo-remanent 
magnetization (trm) relaxation at a given temperature, 
which shows that the short time and long time singular- 
ities are described by different exponents Xk @- It is 
natural to interpret the different levels of the hierarchy 
in terms of length scales, as is actually suggested by the 
RSB solution of the pinned manifold problem [|15j . There- 
fore, one can expect that for any temperature within the 
spin-glass phase, there will be a particular 'critical' length 
scale I for which x — x(£) = 1, that will contribute to 
the correlation and response functions as a function of 
\ogt/ \ogt w . We have reanalyzed some trm data in the 
spirit of the present discussion. We show in Fig. (g) the 
decay of the trm in AgMn at T = 0.75T S @, plotted as 
a function of s(t, t w ) — log(l + t w /t)/ \ogt w , as suggested 
by Eq. (g). This fi gure shows that the rescaling is very 
good at short times, but is violated for large t/t w (cor- 
responding to small s(t,t w )). From Eq. (28), we indeed 



expect to observe the sum of a contribution A4(k < k*) 
from the levels k < k* (that only depends on t/t w ) and a 
contribution from k ~ k* , proportional to s(£, t w ). When 
a = logt/ logt w < 1, the first contribution does not vary 



much, since (t/t w ) 



t 



(a-l){l-x k ) 



is very small for 



large t w , whereas the second contribution is a function 
of a. This suggests that in the short time regime one 
should observe: 



M(t w + t, t w ) = M{k < k*) + m lS (t, t w ), 



(29) 



where mi is a certain prefactor. The quantity ip — 
mi/(M.(k < k*) + mi) measures the relative contribu- 
tion of the levels < k* and around k* in the total decay 
of the signal. The dashed line shown in Fig. ^ is a linear 
fit of the initial decay, as a function of s, from which one 
extracts (in this particular example) tp « 0.2. 
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FIG. 4. Plot of trm data in a AgMn spin-glass at 
T/T g = 0.75, for t w = 300, 1000, 3000 and 10000 seconds. 
The horizontal axis is the variable s(t, t w ) defined in the text. 
The rescaling is very good for s(t, t w ) > 0.3 approximately (or 
logt/ log t w < 0.7), but becomes inadequate for longer times 
t, as more clearly seen in the inset where log s is used. The 
dashed line is an afBne ht of the short time part of the data. 

In all the above formulas, time is implicitcly measured 
in units of a microscopic time To. The value of To is 
not necessarily an individual flip time ~ 10 -12 sec, since 
collective dynamics may exist in the vicinity of the transi- 
tion point (see for a detailed discussion). In Fig. 
||, we have chosen To = 10~ 5 sec. to achieve the best 
rescaling. This value is very close to the one extracted 
from the analysis of Jl^j ■ 

The conclusion is that the TRM data is indeed compat- 
ible with a \ogt/ logt w behaviour for short times. Note 
that this behaviour is tantamount to a logarithmic depen- 
dence of the a.c. susceptibility, or else to 1/ f noise. We 
have insisted in previous papers |7|,0 on the fact that the 
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existence of levels in the vicinity of x = 1 generically leads 
to 1// noise for long times and low frequencies, because 
the contribution of other levels (both above and below 
k*) become negligible in that regime. We show here that 
the degrees of freedom contributing to 1// noise also give 
rise to dynamical ultrametricity. 
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